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1. Introduction 

In the coming era particle physics experiments will probe ever greater energies so nat- 
urally top quark physics will play an increasingly important role in experimental and 
theoretical studies. A major area of study at forthcoming collider experiments will 
be the precision measurement of the top quark mass, which is one of the fundamental 
input parameters of the standard model and gives rise to the leading contributions 
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to its effective potential. In addition, a thorough understanding of top quark physics 
is crucial for the discovery of new heavy particles, most notably the Higgs boson, 
since, due their copious rate of production, top quarks will provide the main source 
of standard model background. 

Although inclusive quantities such as total cross sections are well described by 
fixed-order QCD calculations, experimental analyses require a detailed description 
of the final state. Large logarithmic contributions to differential distributions must 
be resummed to all orders and hadronization effects taken into account. Parton 
shower Monte Carlo simulations include these higher-order corrections by appealing 
to the strongly ordered soft/collinear limit in which higher-order matrix elements are 
represented by universal factors multiplying the lowest order matrix element. This 
leading-log approximation scheme may be recast in a probabilistic form, a Markov 
chain, from which we can attempt to generate events as they occur in nature. 

In this paper we will describe such a prescription for the simulation top quark 
decays, implemented in the new Herwig++ event generator [1,2]. Many issues raised 
here and aspects of the physics we present are also relevant to a discussion of radiation 
from the top quark in its production. In full generality it is not possible to study 
the production and decay phases separately, doing so corresponds to working in the 
narrow width approximation^ i. e. neglecting the width of the top quark. In practice 
this means the width of the top quark, which is of order 1 GeV, should be considered 
infinitesimal and so we should not simulate gluons with energies below that scale. 
This does not present a problem as such a scale is already approximately equal to 
the typical cut-off scale used to terminate the parton shower (the scale at which the 
parton shower hands over to the hadronization model). 

Simulations of this process have been considered in the past, in the older FOR- 
TRAN HERWIG program [3,4]. The simulation which we describe improves on 
this earlier work in a number of ways, partially due to general theoretical develop- 
ments [5,6]. The older HERWIG program was based on a non-covariant showering 
formalism. This meant that, working in the top quark rest frame, the top quark 
could emit no radiation before it decayed. Consequently the entire phase-space had 
to be populated as though all emissions originated from the 6-quark. In addition, 
the older program had grown out of a formalism based on massless emitting partons 
which forced the introduction of an ad hoc angular cut-off on gluon emission, giving 
rise to an unphysical halo of gluons at small angles [7]. 

In this paper we describe an approach based on a new covariant parton showering 
formalism, which naturally includes gluon radiation from the decaying top quark [5]. 
This has the advantage that the region of phase space corresponding to soft gluon 
emission is populated exclusively by the parton shower, instead of the single gluon 
emission matrix element, as in earlier simulations, which required an ad hoc soft cut- 
off [3,4]. Another significant benefit of this new formalism is that it enables a correct 
treatment of the masses of the emitting partons through the use of quasi-collinear 
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splitting functions: no angular cut-off is required. 

Finally we note that the standard coherent parton shower algorithm has two 
important drawbacks. Firstly, because the parton shower generates emissions from 
each leg of the hard scattering (quasi-) independently, each additional emission must 
be uniquely associated to a particular leg of the hard scattering, which can only 
be achieved at the price of having regions of phase space, corresponding to high pt 
gluon emissions, which are unpopulated by the shower. Secondly, the soft /collinear 
approximation to the QCD matrix elements is not a good approximation all over the 
phase-space region populated by the parton shower. Both of these problems may 
be solved by so-called matrix element corrections [8] which ensure that the hardest 
additional gluon emission in the event is distributed according to the exact matrix 
element. We discuss the inclusion of these corrections in the approach of [5] which 
leads to significant improvements in predictions of physical observables. 

In the next section we present the basic parton shower formalism, based on the 
covariant parton shower formalism described in [5]. As we are using the parton 
shower approximation in a range of kinematics and masses that it has not been used 
in before^ our discussion is accordingly detailed. The inclusion of the matrix element 
correction to the decays is considered in section ^ followed by a discussion of the 
results of the simulation in section |. Finally we present our conclusions and plans 
for further developments. 

2. Basic shower formulation for t^bW decays 

Unlike other quarks, the mass and width of the top quark are such that it emits 
radiation and decays before hadronization occurs [9,10]. In this section we describe 
the parton shower approximation to the decay of a top quark to a ly boson and a 
6-quark with additional gluon radiation. Here we are only concerned with the decay 
of the top quark, that is our initial state is an on-shell top quark produced in some 
hard scattering process. 

2.1 Shower variables 

We use the conventions and shower variables described in [5], which we briefly review 
here; these variables, together with an appropriate ordering condition for multiple 
emissions, ensure that the emissions obey the angular ordering of QCD radiation. 

The decay t bW + (n) g is depicted in figure |l[ The parton showers may 
be viewed as cascades of quarks decaying to quark-gluon pairs: ^ qi + ki. The 
momenta of the quarks and gluons in the shower are defined in terms of the following 
Sudakov decomposition 

qi = aip + PiU + q±i, 
ki = qi-i - qi- 

^In the earlier work of [3] only the 6-quark radiated. 
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Figure 1: An example of the decay t — > hW + (n) g. 

As usual in parton shower simulations, prior to simulating gluon emissions the leading 
order decay must be generated. The initial (on-shell) top and bottom quark momenta 
serve to define the basis vectors p needed for each shower. The n and q^i basis 
vectors are defined in the rest frame of the initial top quark with the initial 6-quark 
momentum in the z direction. In this frame the n vectors for top and 6-quark showers 
are chosen to be 

nt = \mt (1, 0, 1) , 
rib = \mt (A, 0, -A) 

where A is given by, 
1 



(2.2) 



^ = —\l'^t + ^tv + "^b - 2m^m2 - 2mlml^ - 2rn?^ml. (2.3) 

t 



The q±i vectors have the form (0, q_Li, 0) in this frame. 

In addition to the Sudakov decomposition ( |2.1J ) we also adopt the following 
variables from [5] for the top quark shower 



~2 _ m^t-Qi ^ _ 



1' (2-4) 

P±i = q±i - Ziq±i-i, pl^ = -pli, 

where p±i is defined to be the relative transverse momentum involved in branching i, 
and qf is the evolution variable. For radiation from the bottom quark the definitions 
in ( p. 41) are unchanged except for the evolution variable, which becomes 



^1 = (2.5) 



These evolution variables qf are defined by close analogy to those of the original 
HERWIG program [4], they correspond closely to the angle between the emitted 
quark and gluon: ordering in qf results in angular ordering of both the top quark 
and bottom quark parton showers. This feature is discussed in more detail in section 
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2.2 Shower phase space 

We begin by considering the decay t — > bW + (n) g, where the gluons are unresolved. 
The partial width for such a decay is given by 

y"dr„ = ^ y d<l>,wd<t>K (27r)^5^ {pt -q,-q^-K) \Mn\^ , (2.6) 

where nit is the mass of the top, pt is its four-momentum and qt/w is the four- 
momentum of the b/W. We have also denoted the matrix element A4n- K denotes 
the sum of the individual gluon momenta ^"^^ ki, while d^K and d^bw are the 
phase-space measures for the gluon and bW systems respectively, 



/ d^w = ]Ii=b,w I 



(2.7) 



The symbol W denotes the region of phase space inside which gluons are unresolvable. 

Assuming the branching picture depicted in figure |l], the full phase space may 
be factorised by repeatedly inserting the identity as integrals over two simple delta 
functions, one such insertion for each gluon vertex: 

/ d<^bwd<^K = nr=i / / d$,^y (27r)^ 6^ (g, - - q^) . . 

/ d$, = / d<i>,A%dQ^ 6 (Q2 - qf) 54 (^^_^ -q,-k,)' ^ 

where go = Pt, Qn = Qb- Given the branching picture in figure |l|, the parton showers 
may be viewed as cascades of two body decays q qg. This interpretation is evident 
from the fact that d$j has the familiar form of a two body phase-space measure for 
a quark of mass Qi and a gluon, albeit with an additional integration. 

Exploiting the Lorentz invariance of the integration measure, we may rewrite the 
two body phase-space integrals for the top quark (p.8D as 



/ 



= I dqfdz, (1 - z,) , (2.9) 



these differ from those of the 6-quark by a factor of Zi due to slightly different 
definitions of the evolution variables (|2.4|J2.5| ). 

2.3 Soft gluon coherence and angular ordering 

The integration limits on qf and Zi in (p.9|) may be inferred from the kinematic 
constraints involved in each splitting, however, as we shall now discuss, dynamical, 
soft gluon interference, effects motivate further restrictions on this phase space and 
hence a modification to the integration bounds implied by the kinematics. 

Using the soft gluon insertion technique [11], and following a similar analysis to 
that in [12], the squared amplitude for a soft gluon dressing the decay t — > bWg, 
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with amplitude A^i, is given by 

/ 



lim \M2\ 



Cf 
+ Cf 
+ Ca 



where 



1 - 



+ 



rii.nj—ni.nf; 



(2.10) 



(2.11) 



' jk 2 (^ifc ^ij) ' 

u is the energy of the soft gluon, rii = Pi/Ei, Uk = k/u, and the sum over colours 
and spins, is understood. The first two terms proportional to Cf are due to (inco- 
herent) radiation from the top quark and bottom quark respectively, the third term, 
proportional to Ca, is due to radiation from the hard gluon (g). Averaging W-j over 
the azimuthal angle, about pi, one readily finds 



2nj.nfc 



rii-Uk 



Vi {rii.nj - UiMk) 



{ni.rij - rii.nkf - ((1 - ni.nuf - vf) 



(2.12) 

with Vi the velocity of particle i. Clearly in the limit Vi,Vj 1, 



where nj = 

we see that the emission from z in a pair of partons ij is restricted to a cone viz 



2ni.nk 



-0 



(2.13) 



where 6i is the angle between the soft gluon and i, and % is the angle between i 
and j. The vanishing of the radiation outside the cone is attributable to destructive 
interference. 

To demonstrate the angular ordering of the radiation from the top quark we 
need to consider the limit in which the directions of the top quark and hard gluon 
approach one another. We also need to consider fj,t>b — > 1 in order to use ( p. 13 )- In 
the dipole rest frame (the rest frame of the W boson) Vt is around 0.7, so the validity 
of working in the limit ft, f — 1 is questionable, we will address this matter shortly. 
Proceeding with these approximations we find = Ut and hence 



1^21 



^ CfWI + CaW!^ + CfWI, + CfWU l-Mi 



t9 



(2.14) 



The first two terms in the radiation pattern (|2.14| ) are due to (incoherent) emis- 
sions from the top quark and hard gluon respectively. The third term in the radiation 
pattern is due to (incoherent) radiation from the &-quark, due to its interaction with 
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t* , by which we mean the amalgamation of t and g (soft gluon emissions from b will 
not 'resolve' the separation oft and g). In the limit Vt, Vb — > 1, the azimuthal average 
of each of the first three terms, of the form W-j, is given by (|2.13|) , so soft gluon 
emissions are restricted to cones whose half angle is equal to that between the two 
particles forming the colour dipole {i and j). 

Similarly, by virtue of ( |2.13| ), the radiation due to the fourth term in ( |2.14| ) is 
restricted to a wide cone lying along the t* direction, reaching out to the 6-quark, but 
with emissions vetoed in the smaller cones, where the incoherent t and g emissions are 
allowed. The radiation due to the fourth term is considered to be the coherent sum 
of emissions from t and g ( p.l4|) and therefore, equivalently, soft, angular ordered. 



wide angle, radiation from the internal, off-shell t* line. This analysis indicates that 
in evolving toward the decay, from the on-shell top quark, the angle of the emissions 
is always increasing. 

We now turn to question the validity of using the Vt,Vb ^ 1 limit in motivating 
this angular ordering. In [7] it was found that equal finite parton masses in {W-j") 
smooth the angular cut-off provided by the step function in ( |2.13|) to the form shown 



in brackets in (p. 121) , also the finite masses screen the collinear singularity in the 



2ni.T\^ prefactor giving rise to the so-called dead- cone effect [7]. 

In the case of the t — > hW decay, assigning velocities to t and h like those in the 
W rest frame, the dead -cone effect is pronounced but the smoothing of the 6 {6th — Gt) 
cut-off does not occur, as one can see from the plot of in figure |. On the other 
hand, in the t bW case, the radiation distribution does not simply vanish outside 
the cone like 9 {6tb — Ot) but rather it can be negative for small 9tb- Therefore, by 
performing conventional angular ordering, we will be making an approximation i.e. 
neglecting this negative contribution. We argue, as in [7], that this is justified given 
that the majority of the radiation is inside the cone: for Vt ~ 0.7, Vb ~ 1, as in figure 
H, we find that the ratio of the integral of (VT/;,) outside the cone to the same integral 
inside the cone is approximately —0.3. 

The same calculation can be performed for the 6-quark in the limit Vb 1 and 
one finds that the emission angles are required to be successively smaller as one 
evolves away from the decay. 

In terms of the qi evolution variable, in the top quark rest frame, in the soft limit 
{1 — z = e), one may show that 

4 = r^-^' (2-15) 
mf 1 — cos Ob 

where 9b is the angle between the soft gluon and the 5-quark (since the top is at 
rest) i. e. starting at q^ = mf and evolving to higher values monotonically increases 
the angle/transverse momentum of the soft gluon emission with respect to the b- 
quark direction. This is consistent with our analysis of the soft radiation pattern 
from the t — ^ bWg decay (|2.14|) , in the W boson rest frame, which is obtained by 
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Figure 2: In this figure we show the azimuthally averaged radiation pattern (Wj*^), for 
vt = 0.65, Vf, = 1. In this plot, the ratio of the integral of outside the cone, to the 

same integral inside the cone is around —0.33. 

a boost in the W boson direction, which will preserve the transverse momentum 
ordering. Therefore implementing angular ordering in the top quark shower simply 
involves bounding gj+i > g^. For the 5-quark the situation is the same as in the 
older HERWIG program namely g^+i < ZiQi. These bounds give rise to nested qf 
phase-space integrals. Finally we note that this ordering in qf is more stringent than 
the ordering implied by the kinematics, one may easily prove that the qf ordering 
is a byproduct of the gj/angular ordering using just the definitions in section ([2T| ). 

2.4 Matrix element approximations 

Thus far we have shown how the phase space for the parton shower may be exactly 
factorized in a Lorentz invariant way. We will now briefly discuss the factorization 
of the n-parton matrix element. 

The principle requirement for factorization of the matrix element is that the 
emissions are either soft or quasi-collinear or both. Such emissions are significantly 
favoured by the underlying dynamics and so the majority are of this type. We aim to 
accurately describe all quasi-collinear emissions and to take into account soft gluon 
emissions through the approximation of angular ordering. 

2.4.1 Emissions from top quarks 

The quasi-collinear limit [13] is defined to be the limit in which and the on-shell 
quark mass squared (w^g) are assumed small compared to n.p but not compared to 
each other. For an n-parton process in which a decaying top quark emits a quasi- 
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collinear gluon, the squared matrix element factorizes according to 

lim \Mnf = jT^^^^^Pm {zn, ql) \Mn-lf , (2.16) 

qn\\k„ [i — Znj 

(2.17) 

where M.n-\ is the matrix element for the process without the additional gluon 
emission. P^g is the quasi-collinear quark-gluon splitting function [13]. 

Another important requirement for our factorization of the matrix element is that 
the gluon emissions are such that the intermediate quark virtualities are strongly or- 
dered. As mentioned previously, the virtualities are naturally ordered by kinematics 
alone. Furthermore, from the Sudakov decomposition ( |2.1| ) and our definitions (|2.2|), 
one can show that imposing angular ordering automatically provides an ordering of 
virtualities which is more restrictive than that dictated by the branching picture in 
figure |l[ For strongly ordered emissions the decomposition ( p.l6| ) may be applied 
recursively, reducing |A^ri|^ to a product of splitting functions that multiply the 
leading-order matrix element. 

Unfortunately, for emissions from the top quark, given our choice of n ( |2.2| ), 
n.p = \rnl. Therefore one may not expect that the quasi-collinear splitting function 
approximates the single emission matrix element very well. Although the choice of n 
is arbitrary, in the top quark rest frame n.'p only depends on the energy components of 
n and, as we shall see shortly (|2.18|) , the quasi-collinear splitting function is invariant 
under rescalings of n. This is may be viewed as a manifestation of the fact that 
collinear enhancements usually take the form of large mass singular logarithms of 
some scale, characteristic of the leading order process, divided by the quark mass, 
whereas in this case the characteristic scale is the quark mass. This is, nevertheless, 
not a problem; the fact that the collinear emissions are suppressed means that only 
soft emissions are enhanced and, as we shall now demonstrate, these are well modelled 
by our approximations. 

From the definition of the quasi-collinear splitting function and the shower vari- 
ables z and g^, we find that the approximation to the quasi-collinear limit of |A^„|^ 
can be rewritten 

lim \Mnf = SnasCpVt^n , (2.18) 

„ 1 / n.ki 2n.qi n.qifnf \ 

Vt,n = - — T — - + ^ - 7 rr^o • 2.19 

2qi_i.ki \n.qi_i n.ki [n.qi_i) [qi^i.ki) J 

In the soft limit the first term in „ is sub-leading and may be neglected (ki is only 
in the numerator), also in the soft limit qi = qi-i, so our approximation to in 
the soft limit is 

^ ''^■qt 
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Neglecting the 6-quark mass n = qb and we may recognise this soft limit Vt^n as 
having precisely the familiar form of the eikonal dipole radiation function, namely, 



T>,n = 4 f A - A 1 . (2.21) 



2 \qb-ki qt-ki 

Plainly, the quasi-collinear splitting function, together with the definition of n and 
the shower variables [z, g^), reproduces exactly the correct soft limit of |7Vln|^ if one 
neglects the mass of the 6-quark. 

Taking into account the 6-quark mass, the differences between the soft approx- 
imation we use ( |2.18| ) and the exact result ( |2.21| ) will become significant if there 



are emissions for which g^./cj becomes small i.e. for emissions near-collinear to the 
&-quark. However such emissions are generated assuming that they were emitted by 
the 6-quark, i.e. they are produced by the 6-quark shower and not the top quark 
shower. 



The high level of agreement between our approximation ( |2.18| ) and the exact soft 



distribution ( |2.21| ) is understandable given the size of the ratio mb/rrit. Despite this 
nice feature, the distribution of radiation which is produced is nevertheless based 
on an approximation, in particular the distribution of high transverse momentum 
emissions is known to be modelled poorly. This matter will be rectified later through 
the soft-matrix element correction, to be discussed in section ^. 

Finally, although, strictly speaking, a higher order effect, earlier studies, neglect- 
ing quark mass effects, have shown that a careful choice of scale for as enables one to 
include parts of the next-to-leading-log contributions [14]. More specifically, for light 
quarks, the soft z ^ 1 limit of the Pqq massless splitting function at the two-loop 
order may be obtained by expanding as (p^j = (1 — zf' g^) about q^ to O (al) and 
rescaling Aj^ . For massive partons we assume that the appropriate argument for 
as is again (1 — z)^ (f, which is equal to p^j + (1 — Zif mf. Although the argument 
of as is no longer the local transverse momentum and although the quasi-collinear 
splitting function contains a term not present in the massless case, in the limit z — > 1 
both the argument of as and the quasi-collinear splitting function tend to their re- 
spective values in the massless analysis and so we can once again expect to reproduce 
the higher order terms. 



2.4.2 Emissions from 6-quarks 

Gluon emissions from 6-quarks are distributed according to the associated quasi- 
collinear splitting function, this is the same as (|2.17| ) but for the replacements 



Zi-m 



2 



rai and qf zqf. As in the case of emissions from the top quark, we may express 
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this factorization in terms of scalar products by using the Sudakov decomposition: 

hm \Mnf = SnasCpVh^n \Mn-if , (2.22) 

Pfen = ^— f^^ + ^-^V (2.23) 
2qi.ki \n.qi_i n.ki q^.ki J ' 

This splitting function accurately reproduces the distribution of gluons emitted at 
small angles to the quark, however, sizeable errors may occur if the approximation 
is used beyond the quasi-collinear limit. 

Evidently emissions close to the n direction are enhanced by the second term in 
( |2.23| ), proportional to n.k^^ . This enhancement is unphysical; as well as being a 
basis vector of the Sudakov decomposition, n is also the (axial) gauge vector used 
in the calculation of the quasi-collinear splitting functions. The manifestation of 
such divergences indicates a breakdown of the quasi-collinear approximation, i.e. 
emissions for which qi-i.k ^ n.ki are beyond the quasi-collinear limit. 

As was shown for the case of the top quark shower, choosing n collinear with the 
colour partner gives a good approximation to the eikonal limit and assigns a physical 
origin to the n.k~^ divergence. Unfortunately, for the case in hand we cannot choose 
n collinear with the top quark since we are in its rest frame and n is light-like. Instead 
we choose n acolinear with qt, thereby maximally separating it from the region of 
phase space into which the shower can emit gluons. In spite of this, results show a 
substantial excess of high transverse momentum emissions from the 6-quark, due to 
the n.k~^ enhancement. This is later rectified by the soft matrix element correction 



procedure described in section 3.1 



Although the soft matrix element correction will compensate any excess emission, 
we would prefer that this was not the default modus operandi, since we only imple- 
ment these process-specific corrections for cases of special interest. By generalizing 
the quasi-collinear splitting function to 

P„-.V„ = P„-C.(f|)^ (2.24) 

and setting n = qt, we can also remedy the surplus emissions and generally improve 
the shower approximation. This modification produces the correct soft and collinear 
limits, it is similar in nature to a soft matrix element correction procedure [8]. Unlike 
the soft matrix element correction, the generalized quasi-collinear splitting function 
Vqq is process independent. In the most basic sense, our splitting function may be 
viewed as a simple merging of the quasi-collinear limit with the full eikonal limit of 
the colour dipole^. However, our generalization is also motivated from more formal 
considerations viz working with a more general class of gauges, which we discuss in 
appendix 0. 

^In this latter respect Vqq differs from the dipole splitting functions of [15]: reproducing the 
eikonal limit of the colour dipole requires a sum of two dipole splitting functions. 
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2.5 Sudakov form factor and shower algorithm 

Proceeding with strong ordering and quasi-collinear factorization as simplifying as- 
sumptions, this corresponds to working in the leading-log approximation. It is well 
known that large logarithms associated to coUinear and soft emissions must vanish 
when the full phase-space is integrated over, this is a consequence of the Block- 
Nordsieck and KLX theorems. We may use this fact to rewrite the integrals over 
the phase space for unresolved giuon emissions U (where we take U to be the region 
z > Zc where Zc is some cut-off) in terms of integrals over the remaining resolved 
gluon phase-space TZ, specifically we use, 

/ dPi {t ^tg)+ I dPi {t ^tg)^0 (2.25) 
Ju Jn 

where 

/ dP, {t ^tg)^ / dq^dz, -^P,, {zi, ql) O {U/U) , (2.26) 

and Q(U/7V) = 1 for gluon emissions in U/Tl and zero otherwise. This gives the 
following expression for the decay width, with n unresolved branchings, 

y dr„ = f[ dPi {t ^ tg)^ ^ J d^,w (27r)' 5^ {qt,^ - q, - qw) | A^o|' ■ 

(2.27) 

Within leading-log accuracy we may neglect the dependence of the final decay, to 
the 6-quark and W boson, on „. Such approximations are a generic feature of all 
leading-log parton shower simulations. With this in mind we may simply rewrite the 
nested integrals over the branching probabilities as 

SdV^^^{-j^dP,{t^tg)Y /dFo. (2.28) 

Summing over all n. the corrections due to unresolved gluon emission exponentiate 
to give, in the leading-log approximation, 



■ u 



exp {- J dPi {t ^ tg)^ J dVo. (2.29) 



Since all soft and collinear logarithmic corrections to the width must vanish on inte- 
grating over the entire phase-space TZ+U, in the leading log approximation the total 
width is simply equal to that given by the leading-order calculation: Fq. Hence, the 
probability that the top quark evolves from scale go to q^^x without emitting any 
resolvable radiation is, 
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the Sudakov form factor. By the same token dPi represents the probability of a 
branching t ^ tg with Zi in the interval [zi,Zi + dzi] and qf in [qi,qi+dqf], for 
resolvable gluons^. The same calculations hold for the 6-quark. 

With these results we find the probability that the top quark evolves from scale 
Qi to q'^ax without emitting any radiation is 

^^'^^^'^"^<^^) - A{qlqf) • ^^-^^^ 

Therefore the probability that there is some resolvable radiation emitted as the top 
quark evolves from scales qf to g^^^ is 1 — 5 {qf, q^ax)^ which may be written as 



2 ""^ 



dg^ 



Qt -"It 



The integrand of (|2.32|) is the probability density that the next resolvable branching 



after q} occurs at scale (f . 

The probability distribution (|2.32|) is sampled using the veto algorithm as de- 
scribed in [6]. Finally we note that, the gluons produced by the top quark and 
bottom quark will in turn produce their own parton showers; these showers are stan- 
dard time-like gluon showers, not subject to any complications, e.g. initial-state 
evolution or parton masses, as is the case for the top quark shower, the details of 
such gluon showering are described in [5]. 

2.6 Phase-space limits - boundary conditions 

We will now discuss the integration limits for the Sudakov form factor starting with 
z. The boundaries on the z variable depend on the (f. For a given (f the z boundaries 
are determined by the requirement that p^^ > 0. In the case of the top quark the re- 
quirement that the intermediate top quark has a mass greater than (mb + niw)^ pro- 
vides a further constraint on z. The implementation of these complicated z bounds 
is straightforward using the veto algorithm [6]. 

The lower bound on (f for both top and bottom quarks is simply if = due 
to the fact that the top quark is initially on-shell and the final 6-quark must also 
be on-shell. The phase-space boundaries in the variable were calculated in [5]. It 
was shown in [5], that for a given qb,max there is an upper bound qt,max, for emissions 
from the top quark, 

{Qt,max - ml) {ql^^^ - ml) = i {ml - m^, + ml + Xf , (2.33) 



^With this interpretation (|2.25|) implies that this same probabiHty dPi is negative in the un- 
resolvable region. Clearly the interpretation of dPi as a probability in this region is not sensible, 
however the minus sign is expected, it is known to originate from the virtual gluon emissions (loop 
contributions), which are of course unresolvable emissions. 
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that partitions the phase-space for the decay t — > bW + g into three regions which 
do not overlap: one region accessible to gluon emissions from the bottom quark, one 
accessible to emissions from the top quark, and a further inaccessible dead region. 

This partitioning is vital to ensure that each phase space point has a unique 
matrix element approximation assigned to it and to avoid double counting of phase- 
space points. As such the three regions do not overlap, this is particularly significant 
for the soft gluon region of phase space, which is most heavily populated. The price 
to be paid for this partitioning is the presence of the dead region, however this region 
is solely comprised of high p_L gluons, so it is sparsely populated, and in any case 
the parton-shower approximation is known to model such emissions badly. 

Given the relation (|2.33| ), between g^^aa, and qlmax^ define two choices of 



phase-space partition, the so-called symmetric choice, 

Qlmax = ml + \ (m? - ml, + ml + \) (2.34) 

and maximal choice 

^Imax = 4 {{mt - mwf - ml) . (2.35) 

In the case of the symmetric choice, the phase-space volume is divided more or 
less evenly between that accessible to emissions from the top and bottom quarks, 
while the maximal choice maximises the volume to be populated by emissions from 
the 6-quark. The two choices of phase-space partitioning can be seen in figure ^. 
Typically we favour the symmetric choice, as the maximal choice involves generating 
high transverse momentum emissions from the 5-quark. 

2.7 Kinematics reconstruction 



Given the boundary conditions and the emission probability distribution (p. 32]) we 



have all we need to generate parton showers in terms of q^^ Zi and 0,. The Sudakov 
variable may be calculated from ( |2.4|) as each Zi is generated, the transverse 
momentum is also calculated at this point. The calculation of the variables is 
more complicated. 

For time-like, final-state, showers we start at the end of the shower, since the 
particles are on-shell there we can simply determine the value of (3i by computing ql 
from ( |2.1|) and setting = m^, where is the on-shell mass squared. Once this 
is done the parent particle's virtuality, and hence its (3i value, follows directly from 
momentum conservation. Applying momentum conservation to each vertex in such 
showers, one can fully reconstruct all of the momenta up to and including that of 
the shower progenitor particle. 

Recall that the gluons radiated by the top quark will produce their own time- 
like showers. The first step in reconstructing the momenta of the top quark as it 
evolves toward its decay, is to reconstruct these showers, i.e. we first reconstruct the 
momenta of the gluons that were radiated by the top quark by applying the procedure 
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Figure 3: In this figure we show the phase space boundaries for the symmetric (left) and 
maximal (right) choices of phase space partitioning, in the Xg,xw plane, where Xg and xw 
are equal to two times the energy fraction of the gluon and W boson in the top quark rest 
frame. The regions Tl and T2 are populated by gluon emissions from the top quark while 
the region labelled B is populated by emissions from the 6-quark. The region labelled D 
is the dead region. For the symmetric choice the phase space volume is divided more or 
less evenly between that accessible to emissions from the top and bottom quarks, while the 
maximal choice maximises the volume to be populated by emissions from the 6-quark. 

outlined above. The momentum of the top quark after each gluon emission can then 
be calculated by momentum conservation, as the initial top quark momentum is 
known (mf,O,0). In this way we completely determine the momenta of the top 
quark and all of the radiation emitted prior to the t — > bW decay. Should the b- 
quark also emit radiation, the resulting jet (s) will be reconstructed in the standard 
way for time-like showers [6]. 

Since the initial b and W boson momenta were generated according to the tree 
level t — ^ bW decay they add up to {rrit, 0,0), rather than the actual momentum 
of the top quark prior to its decay. Furthermore, if the &-quark radiates, forming a 
ft-jet"' the jet reconstruction described above leads to the initial 6-quark momentum 
having a virtuality greater than its mass. 

In order to ensure global momentum conservation we perform a 'momentum 
reshuffling' which smoothly preserves the internal properties of each jet. Let us 
denote the total momentum of the jets radiated by the top quark gisR- The first 
step in the momentum reshuffling involves rescaling the three-momenta of the W 
boson to account for the loss of energy due to gluon emissions from the top quark. 
The second step involves absorbing the component of Qisr transverse to the W, in 
the 6-jet. Finally the momentum of the 6-jet in the direction of the original 6-quark 
is rescaled. This process is sketched in figure ^ and is described by the following 

^In the present discussion we will assume the 6-quark has given rise to a jet since this is the 
general case, the extension to the case where it does not radiate is trivial. 



b b 




Figure 4: In this figure we sketch the 'momentum reshuffling' procedure. Initially the 
top quark decay to a 6-quark and a W boson is simulated, this momentum configuration 
is shown first on the far left. Afterwards some additional radiation is produced from the 
parton shower, represented by the spiral in the central configuration above, clearly this 
configuration does not conserve momentum. The initial 6-jet momenta are rescaled and 
boosted to give the configuration shown on the right hand side: the three-momenta of 
the W boson and 6-jet are scaled down and the momentum of the additional radiation, 
transverse to the VF-boson direction, is absorbed by the &-jet . 



mapping: 




qw = (^/"^^ + P^-p) ^ (A/"^VF + ^2P^-^2P 



'Let + ^?P^ + P?5R±' ^iP " PlSR±) , 



(2.36) 

where p is the initial W boson momentum from the t — ^ bW process, ki and k2 are 
constant rescalings of the 6-jet and W boson momenta. The three vectors pisR\\ and 
PiSR± are the components of gisR parallel and perpendicular to p. Applying energy- 
momentum conservation to the momenta above allows one to determine the value of 
the rescalings ki and ^2. The top quark and the radiation from it are untouched by 
the momentum reshuffling. 

In practice, rather than applying the rescalings and momentum subtractions as 
indicated in ( |2.36| ), once ki and k2 are determined we actually calculate the Lorentz 
boosts which perform a mapping that is exactly equivalent to that in ( 2.36|) . By 



calculating the mapping in terms of these Lorentz boosts, we can conserve momentum 
and preserve the internal structure of the jets by applying the boosts to each particle 
in the 6-quark jet. 

If the top quark does not radiate any gluons before decaying, the three-momenta 
of the 6-jet and W boson system are rescaled such that its invariant mass is equal to 
the top mass. In this special case the problem of momentum reshuffling is identical to 
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that of e+e^ qq events, hence for these cases we use exactly the same momentum 
reshuffling as described in [6]. 

This sequence of boosts and rescalings is designed such that, at least for the case 
of one gluon emission, the 3-body decay kinematics of [5] are reproduced, i.e. should 
the gluon be emitted by the top quark, or the 6-quark, its transverse momentum is 
absorbed by the 6-quark, in the top quark rest frame (as depicted in figure^). 



3. Matrix element corrections 

As stated in section |, the effects of unresolvable gluon emissions have been included 
to all orders through the Sudakov form factor. The master formula and shower 
algorithms generate further resolvable emissions by approximating the full t — > bW + 
(n) g matrix element by a product of quasi-collinear splitting functions multiplying 
the tree level amplitude. Ideally, we wish to include the higher-order effects in the 
master equation as accurately as possible. 

3.1 Soft matrix element corrections 

In the parton shower approximation the probability density that the ith resolvable 
gluon is emitted into [g^, + dg^] , [z,z + dz] is 

dP {z, = ^dz ^^P,, {z, e (7^) . (3.1) 

This approximation works well for the case that the emission lies within the domain 
of the quasi-collinear limit. On the other hand the exact matrix element calculation 
gives us that the probability of a resolved emission is (at least in the leading log 
approximation) 

dP-^- = / dfdz ^/^e (7^) , (3.2) 

J To dzdq^ 

where F is the width of the process t bWg. The differential width for t bWg 
is given in appendix [B.2| . Once again the 'probability' of an unresolved emission can 
therefore be written — dP™ *^, proceeding in the same way as our earlier derivations 
( |2.32|) , we then have the probability density that the ith resolvable gluon is emitted 
into [g^, g^ + dg^], [z, z + dz] is given by the integrand of 

1 d^r / 1 d^r 

^^i^^ exp - / dg dz ^^-r-fr^ 

g2_^ Todzdqt \ Jg2_^ Tq d2;dg^ 



We may generate the distribution in (|3.3|) by simply augmenting the veto algorithm 
that is used to produce (p.32|) with a single additional rejection weight, simply vetoing 



emissions if a random number TZs is such that 

^-^dp 



(3.4) 
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For this to work we require that the parton shower emission probabihty dP al- 
ways overestimates that of the exact matrix element dP^^^■, if necessary this can be 
achieved by simply enhancing the emission probability of the parton shower with a 
constant factor. 

In practice we do not apply this correction to all emissions as most should be 
well approximated by the parton shower. The emissions which are not expected to 
be modelled well are neither soft nor collinear, they have large transverse momentum 
(p_l). In practice we make the ansatz that all emissions which do not have the largest 
p_L of any generated thus far, are considered to be infinitely soft and we only correct 
the set which is comprised of those emissions which have the largest so far, to 
the full matrix element distribution. 

One might be concerned that in this case it is really only proper to apply this 
correction to the final, largest p^ emission, however, in the context of our coherent 
parton branching formalism (angular ordering) the earlier wide angle emission is 
considered too soft to resolve the subsequent, larger pj^ splitting, and is therefore 
effectively distributed assuming that the latter emission did not occur. Under these 
assumptions the correct procedure should be to correct only those emissions which 
are the hardest so far, from distribution ( p.l|) to distribution (|3.3|) by applying veto 



3.2 Hard matrix element corrections 

In addition to correcting the distribution of radiation inside the Tl, T2 and B regions 
populated by the parton shower, we also wish to correct the distribution of radiation 
outside, in the high p^ dead region. We wish to distribute the radiation in the dead 
region according to the full t bWg matrix element i.e. according to 



where b = ml/ ml, w = m^/ml and A is given by (|2.3|) . The full expression for the 
width is lengthy and so we give it in Appendix |B.2| . 

The first step in the algorithm is to generate a point inside the dead region. This 
is non-trivial as the phase-space boundaries shown in figure ^ are rather complicated 
functions of the mass of the 6-quark, W^-boson, z and (f . It turns out that the phase- 
space may be parametrized as functions xw {xg, g^), that is, the Dalitz variable xw = 
2qw.pf/ml may be written as a function ofxg = 2qg.pt/ml and the evolution variable 
q^{to do this one eliminates z in favour of Xg). These functions were calculated using 
the conventions in [5] and are noted in appendix |B.1| . 

The phase-space point is selected by importance sampling the differential dis- 
tribution ( p.lO| ) assuming a {1 + w — xw)~^ behaviour, where a is a parameter 
which may be tuned to improve the sampling efficiency. First we rewrite the integral 



in (U [8]. 



1 d^F 

Fo dxgdxw 




(3.5) 
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over the dead region (D) as 

dxgdxw 



/d (ixgdxw d^g^ = T-oIb dygdyw W {Xg, Xw) 



W{Xg,Xw)= x"{l + w—xw) '^^^ 



dxgdxw^ (3.6) 

= In (1 + w - Xw) , 

and generate yg, yw points assuming they are uniformly distributed between their 
maximum and minimum values in the dead region. Then we keep events with prob- 
ability 

{V < 1) where 

V (yg) = {yg ,max yg,min ) (yw, max yW,min)\y^ i 

(3.8) 



\max 



is the Monte Carlo estimate of the volume of the dead region in yg, yw plane: yg^r. 
and yg^min are the maximum and minimum possible values for yg in the dead region, 
while yw,max and yw,min are the maximum and minimum possible values of yw in the 
dead region, for the generated value of yg. 

Once a pair of values yg, yw has been successfully generated, they may be 
mapped back to the corresponding Xg, xw phase-space point and the momenta of 
the top quark, bottom quark and W boson can be constructed from there. 



4. Results 

In figure ^ we show the Dalitz distributions for the process t bWg as given by the 
parton shower algorithm, including matrix element corrections. Only decays where 
one gluon is emitted are considered, since to do otherwise would constitute a different 
(higher dimensional) phase-space volume to that shown. 

The Dalitz plots show significant clustering of events in the soft Xg ^0 region 
as expected. One can also see a marked clustering of events along the line xw ~ 1-2 
where the energy fraction of the W boson is maximal i. e. where the gluon is emitted 
collinear with the 6-quark. Again, this is to be expected since we know that collinear 
emissions from the 6-quark are enhanced by logarithms of m^/m^. We do not see a 
similar wedge-like clustering of events at low transverse momentum with respect to 
the top quark. Such events would lie along the lower section of the red boundary 
between x^, ~ and Xg ~ 0.75. Clustering along this line would only come from 
collinear enhancements of the form log % where is some hard scale, but in this 

2 

case = mf so there is no such enhancement, only soft enhancements log^ giving 
rise to the observed clustering as Xg —>■ 0. 

In addition to considering the isolated decay process we have also considered 
the process e'^er tt near threshold: ^/s = 360 GeV. By working close to the 
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Figure 5: Dalitz plot for gluonic radiation in top decay. In both plots the soft and hard 
matrix element corrections have been applied, but only one emission has been allowed, a) 
shows the radiation for the symmetric choice of [5] for emission from the top and bottom 
while b) shows the radiation with the scales chosen to give the maximum amount of radiation 
from the bottom quark. The blue (innermost) line gives the limit for radiation from the 
bottom, the green (middle) line from the top and the red (outer) line the boundary of the 
phase space region. 

tt threshold we inhibit the effects of radiation from the production phase of the ti 
pair, thus highlighting the effects of radiation in the ti decays. In analysing these 
events we have worked at the parton level and only considered leptonic W decays. We 
have clustered all final-state quarks and gluons into three jets using the k± clustering 
algorithm [16], taking care to omit the W decay products from the clustering. Events 
for which the minimum jet separation is less than AR = 0.7 (Ai?^ = At]'^ + A0^), as 
well as events containing a jet with transverse energy less than 10 GeV, are excluded 
from the analysis. For these events we have plotted the distributions of the jet 
separation AR and the logarithm of the jet measure y^, where 

1/3 = ^ min,, (min {E^ E^) (1 - cos %)) (4.1) 

is the value of the jet resolution parameter for which the three jet event would be 
seen as a two jet event. This analysis is the same as that performed in two earlier 
related works [3,17]. 

In figure ^ we show differential distributions with respect to the jet separation 
of the closest pair of jets in the event (Ai?), for the cases of the symmetric and 
maximal phase space partitions ([2.34[ , p.35| ) . Figure [| shows that the matrix element 



corrections have significant consequences for the AR distributions. We see that when 
the soft matrix element correction is applied, the distribution is more peaked for small 
AR and tails off more quickly than the distributions obtained without it. This is 
to be expected, it is a softening of the distribution, it is indicative of the fact that 
the soft matrix element correction is vetoing a number of high Pt emissions, not well 
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Figure 6: where /S.E? = /S.rf + Ac/)^ for three jet e+e" tt events at -y/s = 

360GeVwith and without matrix element corrections. On the left we show the distributions 
obtained for the case that the phase-space volume populated by the parton showers from the 
top and bottom quarks is almost the same size (the symmetric phase space partition ( 2. 3-^1 )), 
while on the right we show the distributions obtained for the case that the shower from 
the 6-quark populates most of the phase space choice (the maximal phase space partition 
( 2.35 )). In each plot the black line corresponds to the parton shower approximation, the red 
line corresponds to the parton shower including hard and soft matrix element corrections, 
while the green/blue lines respectively correspond to including only the hard/soft part of 
the matrix element corrections. The magenta line is obtained using the standalone parton 
shower but with Pqq replaced by the generalized quasi-collinear splitting function Vqq (2.24). 



modelled by the standalone parton shower: such hard emissions naturally give rise 
to more widely separated jets (larger Ai?). 

The degree to which the soft matrix element correction affects the Ai? distribu- 
tion suggests that the standalone parton shower does not model the number of these 
high pt emissions well: similar distributions were observed in version 5.9 of the older 
HERWIG program and were since considered to be the result of a bug (showering 
in the wrong reference frame). In our case the effect shown is understood to be a 
genuine artefact of the covariant parton shower formalism and may be traced back 
to the form of the quasi-collinear splitting function for 6-quark emissions. In section 



2.4 we noted that, as the angle between the reference vector (n) and the emitted 



gluon {k) decreases, emissions will be enhanced in the n direction (particularly soft 
emissions). This spurious enhancement results from choosing n not equal to precisely 
the momentum of the colour partner of the 6-quark. This also explains why the dif- 
ferences between the results with and without matrix element corrections are more 
pronounced for the maximal phase space partition than the symmetric partition. On 
the contrary our results including the soft matrix element corrections compare well 
with those obtained in the earlier FORTRAN HERWIG program. 



In section 2.4 we also proposed a new generalized quasi-collinear splitting func- 
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Figure 7: idi^^ where 2/3 = 1 i^mj (min (^Ef, sfj (1 - cos%)) for three jet e+e ^ 
tt events at = 360GeVwith and without matrix element corrections. As in figure |6|, on 
the left we show the distributions obtained for the symmetric phase space partition ( |2.34D , 
while on the right we show the distributions obtained for the maximal phase space partition 
(2.35. The histograms are coloured in the same way as for figure H. 



tion Vqq (see also appendix which improves on the usual quasi-collinear splitting 
function by reducing to the eikonal dipole function in the soft limit. From a practical 
point of view the introduction of this splitting function is akin to a soft matrix ele- 
ment correction. The Ai? distributions support this assertion; the standalone parton 
shower is greatly improved by the use of this new splitting function. 

As with figure ||, in figure we show differential distributions with respect to the 
jet measure y^, for the cases of the symmetric and maximal phase space partitions, 
for all possible combinations of matrix element corrections. We also show the same 
distribution obtained using the generalized quasi-collinear splitting function in the 
basic parton shower. When the soft matrix element correction is applied, the distri- 
bution shifts toward smaller y^ values since some high emissions will be vetoed by 
it. Conversely the hard matrix element correction supplies more events with high pt 
emissions and so the distribution shifts toward higher y^, values, since these events 
will not require such fine resolution to distinguish three jets from two jets. As with 
the Ai? distributions, the introduction of the generalized quasi-collinear splitting 
functions, V^g, gives a substantial improvement of the basic parton shower, softening 
the jet structure. 

In figures || and ^ we see the effects of varying the choice of phase space parti- 
tioning, as well as the gluon mass parameter, on the Ai? and y^ distributions. The 
distributions are less sensitive to changes in this parameter when the matrix ele- 
ment corrections are applied, compared to the case of the standalone parton shower. 
This is to be expected given that the quasi-collinear splitting functions {Pqq) are 
known to produce an excess of emissions as one approaches the acolinear direction. 
This explains the variation of the distributions on shifting the phase space bound- 
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Figure 8: Here we show the distribution as in figure ^. On the left we show the 

distribution obtained from the parton shower approximation and on the right we show the 
same distribution including all matrix element corrections. In each case we have shown 
how doubling the gluon mass and/or varying the choice of phase space partition affects the 
result. 
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Figure 9: Here we show the ^37^5^^ distribution as in figure 0. On the left we show 
the distribution obtained from the parton shower approximation and on the right we show 
the same distribution including all matrix element corrections. In each case we have shown 
how doubling the gluon mass and/or varying the choice of phase space partition affects the 
result. 

ary, moreover the majority of these spurious emissions will be soft, resulting in a 
heightened sensitivity to the gluon mass parameter. 

5. Conclusions 

At the beginning of this paper we presented a theoretical framework for the sim- 
ulation of QCD radiation emitted in top quark decays. In section || we described 
the basic parton shower formulation. Although this formalism is the product of al- 
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most thirty years of evolution, it has not previously been applied to the extremes of 
kinematics and mass scales found in top decay. In doing so we find that the latest 
covariant parton shower formalism of [5] works well under these conditions and we 
have gained new insight regarding its use. In particular, we introduced a general- 
ization of the quasi-collinear limit [13, 15] for which the q ^ qg splitting function 
correctly reproduces the eikonal limit. 

In section |^ we have described the soft matrix element correction which fixes the 
distribution of the hardest emission to be that of the leading-order matrix element 
for the decay t bWg, in a manner respecting colour coherence. We also describe 
the hard matrix-element correction which simply populates the dead region of phase 
space (not populated by the parton shower) according to the same single gluon 
emission matrix element. 

The results of our program show good agreement with those of the FORTRAN 
HERWIG program for the AR and 1/3 observables, provided that either the soft 
matrix element correction is applied [3], or our proposed generalization of the q qg 
quasi-collinear splitting function is used (|2.24| ). The associated Dalitz plots obtained 
for isolated t bWg decays also meet with our expectations. Distributions obtained 
using the standalone parton shower with the conventional quasi-collinear splitting 
functions, suffer from an excess of high pt emissions from the 6-quark. This excess 
in the original basic shower formulation is well understood to arise from the choice of 
reference vector (n) required for the ^-quark shower. Future versions of the Herwig++ 
program will use the new splitting function by default. 

The simulation we have presented improves on that of the older HERWIG pro- 
gram in a number of ways through the use of the new covariant parton shower 
formalism. The new formalism brings with it the use of the quasi-collinear Altarelli- 
Parisi splitting functions, which lead to a natural screening of collinear singularities, 
allowing us to dispense with the ad hoc angular cut-off, which was responsible for 
the spurious dead- cone halo. Moreover, in this paper we propose a generalization of 
these splitting functions, leading to further improvements in the modelling of soft ra- 
diation. We have also been able to generate radiation from the top quark in its rest 
frame, thereby populating infrared regions, corresponding to soft gluon emissions 
from the top quark, with the parton shower. Previously such regions were popu- 
lated according to a single emission matrix element correction with an arbitrary soft 
cut-off [3]. 

It is clear from our discussion that this work may be readily extended to decays 
of other heavy particles, in particular squarks and gluinos. Should supersymmetry 
be realised in nature, squark and gluino decays will give rise to a significant amount 
of activity in the LHC experiments, which we will need to simulate. Moreover, as 
with the top quark today, measuring the masses of these particles will be a major 
area of study, requiring accurate simulations of their decays. It is our intention to 
make this extension in a future version of the Herwig++ program. 
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A. Generalizing the q ^ qg quasi-collinear splitting function 

The quasi-collinear limit is that in which q± becomes O (mg) and small [15]. This 
region can be identified by the uniform rescalings q± Xq±, nig XrUg and exam- 
ining the limit A ^ 0. In the case of a quasi-collinear gluon emission from a quark, 
in an arbitrary n particle process, the squared matrix element factorizes as 

lim \Mnf = , P,, \Mn-lf , (A.l) 

where z = n.Qi/n.qi^i. 

Note that, as well defining the Sudakov decomposition, the reference vector n 
also specifies a choice of axial gauge. Restricting the form of n imposes restrictions 
on the choice of axial gauge. To make calculations easier, n is typically chosen to be 
lightlike and perpendicular to k±, as in [13]. The resulting splitting kernel (|A.2|) is 
invariant under gauge transformations respecting these constraints {n^ = n.k± = 0) 
i.e. Lorentz boosts of n, in the n direction. Invariance under more general Lorentz 
transformations requires k± = 0; in this sense kj_ may be regarded as a measure of 
the gauge dependence of approximation ([A.l[ ). 



Parton shower simulations use approximation ([A.l[) beyond the collinear limit, so 



strictly speaking the approach is not gauge invariant. This is not a problem provided 
that gauge dependent contributions are sub-leading. However, in the case of the b- 
quark shower we have noted an excess of emissions at the edge of the shower phase 
space, due to an unphysical singularity, proportional to n.k~^, where n is collinear 



with the W boson. Furthermore, as discussed in section [2]J, in the limit of soft 
emissions the splitting kernel ([A.l|) does not reduce to the correct eikonal dipole 
radiation function. 

What we require is that the splitting kernel reproduces the correct collinear and 
(ideally) soft limits of the matrix element, without introducing any other singular 
terms. Since the behaviour of matrix elements in these limits is universal, any split- 
ting kernel satisfying these criteria will be, at least to leading order, gauge invariant. 

It turns out that all of these problems can be solved by simply relaxing the 
restrictions on the gauge vector. It is well known that the eikonal limit of a colour 
dipole can be calculated by considering gluon emission from just one quark, provided 
that the gauge vector is equal to the momentum of the colour partner. In order to 
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reproduce the correct soft behaviour we should therefore always set the gauge vector 
equal to the four-momentum of the colour partner of the emitter. Although one can 
often obtain a good approximation to the eikonal limit by choosing n to be lightlike 
and collinear to the colour partner, if the colour partner is heavy, as in top decay, 
this approximation fails. 

In any case, this prompts us to consider the case that n is massive. If we do this 
we find that the steps leading to ([A.l| ) now give, up to sub-leading terms, 

with z defined as z = n.qi/n.qi_i. The additional term in the generalized quasi- 
collinear splitting function (Vgg) directly arises from the term in the gluon polar- 
ization sum. The new splitting function Vgq ( |A.3| ) reproduces the correct soft and 



collinear limits in the shower phase space, it contains no unphysical divergent terms. 

From a purely pragmatic point of view Vgq can be considered as a kind of global 
soft matrix element correction and it is seen to have similar effects on physical ob- 
servables (see section Crucially, unlike the soft matrix element corrections, the 
generalized quasi-collinear splitting function is process-independent. 

B. t—^bWg phase-space and matrix element 

In this appendix we give the matrix element for t bWg decay and the phase-space 
parametrization, these are necessary for the hard matrix element correction discussed 
in section ^. Both the phase-space and matrix element are parametrized in terms of 
the Dalitz variables 

Xi = —, (B.l) 



which, in the top quark rest frame, are equal to two times the fraction of the top 
quark's energy carried by particle i. We also define the following ratios of masses for 
convenience: 

b = , w = — (B.2) 

B.l Phase space 

The Dalitz variables x^, were calculated in [5], in terms of z and the top quark 
evolution variable kt = q'^/mf, assuming the gluon was emitted by the top quark, as 
being 



/ ~ \ _ l+w-b-{l-z)i^t-^/(l+w-b-{l-z)Ktf-Aw{l-{l-z){Kt-l))z 
Xw [Z, K,t) — 2z 

l+w-b-{l-z)iit+y/il+w-b-{l-z)iitf-'iw{l-{l-z){K.t-l))z m 3) 

+ 2(l-(l-z)(St-l)) 

Xg {Z, Kt) = {1- z) Kt, 
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where b = ml /ml and w = m^r/ml. We may completely eliminate z from Xw to 
give 

Xw {xg, kt) = 2(Kt-xg) ifiti^ + w- b-Xg)-A (xg, kt)) 

+ 2(k^+xlil-!i,)) {l+W-b-Xg)+A {Xg, kt)) , 



K+j [Kt - /«- 

2 



A {Xg, kt) — ^ {Xg — Xg^^ [x g " X g S) {^t 

Xg^ = 1 - (^y/w ± Vb 

k± = 2Xg (xg ± ^(1- ^) {Xg-Xg^) {x g - X g 
(l-«l)(l±V^)-fe(lT0«) 



(B.4) 



^9± l~w 



The expression for xw given in ([B.4| ) enables one to draw lines of constant kt in the 
Xw, Xg plane. 

Repeating the procedure for the case that the gluon is assumed to originate from 
the 6-quark gives 



^g± 



{xw, Kb) = '2-Xw - z± [xw, kb) a/x^ - 4w 

-li'^+iT^) {2-xw-V^^^), (B.5) 
z±{xw,kb) = ^ (^kb± ^/kl - 4:kb{l + w - b - Xw) 



Inverting (p.5| ) to obtain xw as a function of Xg involves a high-order polynomial 
requiring a numerical solution, neglecting the 6-quark mass an analytic solution 
becomes possible. 

B.2 Matrix element 

The matrix element for the decay t bWg was given in [3, 5] assuming a massless 
&-quark. We calculate the squared matrix element, without neglecting the 6-quark 
mass and find 

1 d^r _ asCp 1 
To dxgdxw 1" XxwXg 

bx^ 

■J^ + {l-W + b)Xg-Xw{l-Xg) -Xl ^^^^ 

+ l+n.-2j-b{2~w~b) (1 + 2W + b) {xw - Xg)^ + 2wXwXg)^ 

Xw = l + w — b — Xw 

where again, b = mf/mf and w = m^/ml. Setting 6 = 0, our expression ([B.(S| ) 
easily reduces to those given in [3,5]. 
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